Composite-fermion crystallites in quantum dots 
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The correlations in the ground state of interacting electrons in a two-dimensional quantum dot 
in a high magnetic field are known to undergo a qualitative change from liquid-like to crystal- 
like as the total angular momentum becomes large. We show that the composite-fermion theory 
provides an excellent account of the states in both regimes. The quantum mechanical formation 
of composite fermions with a large number of attached vortices automatically generates omposite 
fermion crystallites in finite quantum dots. 

PACS numbers: 71.10.Pm,73.43.-f 



The system of interacting electrons confined to a two 
dimensional quantum dot and exposed to a strong mag- 
netic field has been a subject of intense theoretical st udy 
for over two decades 0,11 SSHBIlHSGilllll 
, 0| , both because such quantum dots have been 
realized and studied in the laboratory 0, 0, , and 
because of the possible relevance of its physics to the 
fractional quantum Hall effect (FQHE) |l8| ■ For systems 
with small numbers of particles, exact diagonalization 
can be performed for parabolic quantum dots in the limit 
when the cyclotron energy is large compared to the con- 
finement potential, which shows that the ground state 
energy as a function of the angular momentum (L) has 
a rather rich structure. In particular, downward cusps 
appear at certain values of L, which are consequently es- 
pecially favorable. An understanding of the correlations 
in the quantum dot state that underlie this physics is 
one of the central questions for this system. One would 
also like to know how this ties into our understanding of 
the FQHE, obtained in the thermodynamic limit without 
confinement. 

The approach based on the formation of composite 
fermions has been demonstrated to be successful in a 
range of L values 0, E| . Specifically, a mean-field type 
description, in which the composite fermions are taken as 
non-interacting particles at an effective angular momen- 
tum L* , with their mass or the cyclotron energy treated 
as a phenomenological parameter, predicts cusps in the 
energy at certain magic angular momenta, which are in 
agreement with the actual cusp positions in exact diago- 
nalization studies ||. However, discrepancies appear at 
large L [lfl H^ : here, the actual cusps occur at regu- 
lar intervals in L, which has been interpreted in terms 
of classical crystal-like states [l^. At large L, the re- 
pulsive interaction thus appears to stabilize a (rotating) 
crystal rather than a liquid, which is believed to signal 
a breakdown of the composite-fermion (CF) description. 
That would not be unexpected, because large angular 
momenta correspond to small filling factors (to the extent 
the filling factor is a meaningful quantity in a finite quan- 
tum dot), and it is known, for infinite two-dimensional 
systems, that the CF liquid gives way to a Wigner crystal 



at sufficiently small fillings. 

To understand the nature of the breakdown of the com- 
posite fermion description, a matter of great interest also 
in the context of the FQHE, we have undertaken an ex- 
tensive study of finite systems at large angular momenta. 
Our investigations, however, have led to a surprising con- 
clusion: Even though the naive mean-field interpretation 
in terms of free composite fermions becomes invalid at 
large L, the microscopic composite fermion theory, de- 
fined through wave functions, continues to give a very 
good description down to the largest L studied to date. It 
provides an accurate approximation for the ground state 
wave function and the ground state energy at every single 
L in the wide range studied, and correctly reproduces all 
cusps in plot of the ground state energy vs. L. Taken to- 
gether, these results constitute a detailed verification for 
the validity of the composite fermion theory for quantum 
dots even at very low fillings. 

This may appear to be at odds with the classical 
crystal-like correlations found in exact diagonalization 
studies 0, While both composite fermions and 

the crystal are generated by the repulsive interaction be- 
tween electrons, the feeling has been that one excludes 
the other. A notable aspect of our work is the unexpected 
finding that not only is there no logical inconsistency be- 
tween the simultaneous formations of composite fermions 
and crystal-like structures at low fillings, but that the 
state is very well described as a crystallite of composite 
fermions. 

The Hamiltonian of interest is 

^ 2m b V 3 c 3 J 2 3 ^ er 1k V ; 

j 3 j<k 

where mj, is the band mass of the electron, ujq is a mea- 
sure of the strength of the confinement, e is the dielectric 
constant of the host semiconductor, and rjk = \rj — r&|. 
We will specialize to the case of very large magnetic fields, 
when uj c = eB/mbC » ujq. Only the lowest Landau level 
(LL) is relevant in this limit. In that limit, at each an- 
gular momentum the eigenenergy neatly separates into a 
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FIG. 1: The exact interaction energy V (dashes) is given as a 
function of the angular momentum L for TV = 6 particles. The 
dots are the predictions of the CF theory. Different panels 
correspond to L regions where composite fermions of different 
flavors are relevant. (In each panel, for L related to L* — —6, 
the lower order CF gives a smaller D* .) The energies are 
quoted in units of e 2 /el. 



confinement part and an interaction part: 
E(L) = E C (L) + V(L) 



(2) 



where E C (L) = (Ti/2)[Q, — uj c }L, relative to the lowest 
LL, with £l 2 = lo 2 c + 4wq, and V(L) is the interaction 
energy of electrons without confinement, but with the 
magnetic length replaced by an effective magnetic length 
given by £ = \fW[rn^i. In the following, we will only 



consider V as a function of the angular momentum L; 
it must be remembered, however, that the confinement 
part must eventually be added to determine the global 
ground state. The exact ground state in each L sector is 
obtained by either matrix diagonalization or the Lanczos 
method. The largest system we have studied has a Fock 
space dimension of 509,267. 

In the CF theory |3l|g,|l9j, the interacting state of elec- 
trons in the lowest LL at angular momentum L is mapped 
into the non-interacting electron state at L* , where 

L = L* + pN(N — 1) , (3) 

N is the number of electrons, and p is an integer. The 
wave functions are related as 

*«=7>n^- z *) 2p <i>f ■ (4) 

j<k 

Here $^ are the wave functions for non-interacting 
electrons at L* (which in general occupy several Lan- 
dau levels), a = 1, 2, • • • , D* labels the different states, 
Zj = Xj — iyj denotes the position of the jth electron, 
V indicates projection into the lowest LL, $„ are basis 
functions for interacting electrons at L, and D* is the 
dimension of the CF basis. We will restrict $ Q to states 
with the lowest kinetic energy at L* , and choose p so as 
to have the smallest dimension for the basis. The com- 
posite fermions carrying 2p vortices are labeled 2p CF's. 
At certain values of L, the above prescription produces 
only one state (D* = 1), which is the CF theory's an- 
swer for the ground state. In the notation of Ref. ||, 
this is a compact state, denoted by (No, Ni, ■ • ■), with 
Nj composite fermions compactly occupying the inner- 
most angular momentum orbitals of the jth CF level. At 
other values of L, when there are many CF basis states 
(D* > 1), we diagonalize the Coulomb Hamiltonian in 
the CF basis to obtain the ground state, using methods 
described earlier 0,H3. For any N, there are many val- 
ues of L where the CF theory gives a unique answer, but 
in general, D* increases with N. 

We have carried out exact diagonalization for N < 8. 
It is possible to to go to arbitrarily large values of L 
within the CF theory, but available computer mem- 
ory has restricted our exact diagonalization study to 
L < 135, 117 and 111 for N — 6, 7 and 8, respectively. 
Figure ^ shows the comparison between the CF theory 
and the exact energy as a function of L for N = 6 par- 
ticles, demonstrating that the CF theory predicts every 
energy accurately and obtains all cusps faithfully [21|. 
The comparisons for other values of N are similar. 

We next come to the connection between our approach 
and the crystalline correlations found in previous studies. 
Such correlations are not manifest in the density, which 
is rotationally symmetric for the ground state with a def- 
inite angular momentum, so one must compute the pair 
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FIG. 2: The upper plot shows the pair correlation functions 
for the exact ground state for TV = 6 particles at L = 95, and 
the lower plot displays the prediction of the CF theory. The 
exact ground state is a linear superposition of 69,624 basis 
states, whereas the 6 CF ground state is given by a unique 
wave function denoted by (4,2). The "missing peak" on the 
outer ring indicates the location R of the fixed particle. 



correlation function 



fl (r)~ J f[d 2 rj \^(r,R,r 3 ,---,r N )\ 

i=3 



(5) 



which is the probability of finding a particle at r while 
holding one particle fixed at R. We first compute the 
density as a function of the distance from the center, and 
fix R somewhere on the maximum density ring. As an 
example, Fig. |2 (upper panel) shows the exact pair cor- 
relation function for 6 particles at L = 95, which clearly 
illustrates crystallite formation. The CF theory gives a 
unique state here, (4,2), which has the wave function 
(with N = 6 and 2p = 6) 



# = VA 



N 



i=3 



i-3 



no%--^) 2p e~ E <- iz!|2/4 

(6) 

where A denotes an antisymmetric Slater determinant. 
The projection is accomplished by the method outlined 




FIG. 3: Pair correlation functions for the (17,2) configura- 
tion for N = 19 composite fermions of three different flavors, 
namely 2 CF's, 10 CF's, and 20 CF's (from top to bottom). The 
missing peak on the outer ring indicates the location R of 
the fixed particle. The three figures are plotted on the same 
length scale to illustrate how the size grows with L. 



in the literature [IJ. The energy of this wave function 
is 1.95535(15)e 2 /e£, which compares well with the exact 
energy 1.95061e 2 /e£ The overlap with the exact ground 
state is 0.902. The pair correlation function for this wave 
function is also shown in Fig. [21 (lower plot). It is re- 
markable that the single wave function from the CF the- 
ory provides a good qualitative and quantitative descrip- 
tion of the actual ground state wave function in Fig. |3 
which is a linear combination of 69,624 basis functions. 
We have studied other "cusp states" and found similarly 
good agreement. 

The composite fcrmion wave function at any angular 
momentum L is obtained from non-interacting electron 
wave function at an effective angular momentum L* by 
multiplication by an appropriate Jastrow factor, in ex- 
actly the same way as wave functions are written for bulk 
fractional quantum Hall states. It is crucial to note that 
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no physics of any kind of crystal is put into these wave 
functions by hand at any stage. For example, referring 
to Fig. |21 the parent electron wave function at L* = 5, 
namely A[- ■ •] in Eq. ©, has no crystalline correlations in 
it. It describes a simple non-interacting system with sin- 
gle particle orbitals {0, 0}, {0, 1}, {0, 2}, {0, 3}, {1, -1}, 
and {1,0} occupied, where {n,m} denotes the angular 
momentum m orbital in the nth LL. (The allowed values 
are n = 0, 1, • • •, and m = — n, —n + 1, • • •.) However, 
when the uncorrelated wave function at L* is multiplied 
by an appropriate Jastrow factor to convert electrons into 
composite fermions, the resulting wave function describes 
a correlated electron liquid for a certain range of L values, 
but for sufficiently large angular momenta, it produces a 
state with crystalline correlation; in other words, the for- 
mation of composite fermions automatically causes the 
formation of a crystallite at large L. It is stressed that 
the same class of wave functions describes both the liq- 
uid and the crystal in quantum dots in different parame- 
ter regimes. The emergence of such classical crystal-like 
structures in the intrinsically quantum mechanical state 
of composite fermions is striking. 

To investigate the evolution with the number of at- 
tached vortices, we show in Fig. [3] the pair correlation 
function for N — 19 composite fermions in the state 
(17,2) (with wave function given in Eq. EJ) for three fla- 
vors of composite fermions, with 2p = 2, 10, and 20. 
The crystal becomes better defined with increasing 2p. 
In this case, we do not have exact results to compare 
with (the dimension of the Fock space is astronomically 
large; even for the smallest system in Fig. namely 19 
particles at L — 477, it is estimated to be 10 14 ), but the 
structure at large L is consistent with classical consider- 
ations that predict a three-ring structure with 1, 6 and 
12 electrons located on inner, middle and outer rings, 
respectively j^. Interestingly, when we hold a particle 
in the middle ring fixed, the hexagonal correlations of 
the middle ring become prominent while the outer ring 
becomes more or less uniform. That is consistent with 
earlier studies that show that the radial correlations be- 
tween different rings are weaker than the angular correla- 
tions within a ring, producing, for example, two different 
melting transitions |23j. 

A recent article by Yannouleas and Landman ^3] has 
asserted that the CF theory fails to produce the cusp po- 
sitions for large L, especially beyond L > 75 (for N = 6). 
That assertion is not borne out by our results, which 
demonstrate that the CF theory, when taken beyond the 
most naive mean field picture, predicts each cusp cor- 
rectl y a t least for angular momenta up to L = 135. 
Ref. [l2j has studied an alternative trial wave function 
based on an analogy to the classical crystal ground state 
in a quantum dot. A thorough comparison of the CF the- 
ory and the approach of Ref. is not possible because 
the latter obtains wave functions and energies only for 
certain special values of L. For N = 6, Ref. [l2j explicitly 



quotes energies from their approach for seven values of 
L in the range 75 < L < 135. For these seven angular 
momenta, the CF theory gives lower energy in every case 
except at L — 135. In future, it would be interesting to 
compare the two methods for a larger range of L and N 
to ascertain their respective regimes of validity. 

One may ask what implications our results have for the 
nature of the bulk FQHE state at small fillings. While the 
finiteness of the systems does not allow us to draw a firm 
conclusion regarding whether the thermodynamic state 
at some filling is a liquid or a crystal, the results do indi- 
cate that the CF state possesses a substantial short-range 
crystalline order at small v. This raises the interesting 
issue of whether the actual Wigner crystal at low fillings 
is a crystal of electrons or of composite fermions [24| . 

The composite fermion description of correlated elec- 
tron states in quantum dots at high magnetic fields has 
several appealing features: It gives a unified theory, ap- 
plicable to angular momenta spanning both the liquid and 
the crystal phases; it produces microscopic wave func- 
tions also for those angular momenta where no classically 
stable crystal is available; it obtains the correct crystal 
shape and the ring structure without an explicit consid- 
eration of the classical solution; and finally, through CF 
theory, the understanding of the quantum dot physics 
dovetails nicely with that of the FQHE. Before ending, 
we note that the accuracy of the CF wave functions can 
be improved straightforwardly and systematically by en- 
larging the basis at L* slightly [25J; the goal in this work 
was to demonstrate that the zeroth order theory itself 
works very well. 
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